Stochastic pumping of heat: Approaching the Carnot efficiency 
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Random noise can generate a unidirectional heat current across asymmetric nano objects in 
the absence (or against) a temperature gradient. We present a minimal model for a molecular- 
level stochastic heat pump that may operate arbitrarily close to the Carnot efficiency. The model 
consists a fluctuating molecular unit coupled to two solids characterized by distinct phonon spectral 
properties. Heat pumping persists for a broad range of system and bath parameters. Furthermore, 
by filtering the reservoirs' phonons the pump efficiency can approach the Carnot limit. 

PACS numbers: 05.40.-a 63.22.-m, 44.10.-K 05.70.Ln 



Heat flows spontaneously from objects at high tem- 
perature to objects of low temperature. Directing heat 
against a temperature gradient requires the application 
of an external perturbation or the design of a spe- 
cial system-bath initial condition. A generic nano-scale 
heat pump model includes a molecular element bridging 
two solids. An external force modulates the subsystem 
(molecular) energetics, leading to a net exchange of heat 
between the baths against, or in the absence, of a tem- 
perature gradient. A spatial asymmetry should be built 
into the system in order to define a preferential direction- 
ality. This process is closely related to the ratchet effect 
where a particle current is catalyzed in the absence of 
a dc voltage drop due to the asymmetry in the spatial 
potential and the influence of time dependent forces [l[ . 

Recent studies have carefully analyzed the operation 
principles of prototype classical @, y|, [1] and quantum 
H, H, 0, H, M, E3, El thermal machines, seeking to 
optimize performance. Typically, in these schemes a 
carefully- shaped external force periodically modulates the 
levels of the nano object, leading to the pumping op- 
eration. In electronic systems a Brownian electron re- 
frigerator based on selective tunneling through a metal- 
insulator-superconductor junction was recently proposed 
Fluctuating electric fields can also do chemical 
work, driving ions against an electrochemical potential 
difference [131 ] . 

In this paper we demonstrate facile pumping of heat 
in the form of vibrational energy in response to random 
noise. The following generic effect is discussed: Heat can 
be pumped between distinct solids when the bridging ob- 
ject suffers stochastic fluctuations of its energies, given 
that the spectral density of at least one solid strongly 
varies within the noise spectral window. Furthermore, 
using reservoirs with suitable energy filters, (Einstein 
solids) one can construct a heat pump that works close 
to the Carnot limit. Thus, we demonstrate that not only 
can one utilize random noise and pump heat against a 
bias ("Maxwell demon"), the efficiency of this machine 



may approach the maximal value. 

Our minimal model consists a Kubo oscillator [3] 
bilinear ly coupled to two macroscopic reservoirs. The 
model could be realized by adsorbing a molecule on 
a dielectric surface, blocking charge transfer, while al- 
lowing for heat exchange between the modules. The 
total Hamiltonian includes the following contributions, 
H = H (t) + H B + V L + V R . H includes the isolated 
subsystem, a Kubo oscillator, 

H (t)= E n (t)\n)(n\; E n (t)=[E^+e n (t)].(l) 

n— 0,1. .oo 

E^ = uBq is the static energy of level \n) , e„(i) is a time 
dependent stochastic contribution. We assume that the 
average over stationary fluctuations vanishes, (e(i)) e =0, 
while higher moments survive. This local oscillator is 
coupled to two thermal baths of temperatures T v = 1 //?„ 
(ks = 1), represented by sets of independent harmonic 
oscillators (h = 1) 



H B = ^2uj k bl k b^ k - 



(2) 



v.k 



b\, k , b V: k are the bosonic creation and annihilation oper- 
ators respectively for the mode k of bath v = L,R. The 
interaction between the subsystem and the bath is taken 
to be bilinear (displacement-displacement) [l5| . 



V„ = F v ^2 Cn\n)(n — l\+c.c. 

n— l,2..oo 

F v = Y, X ^ b lk + b,,k), (3) 



where c n = y/n, and the system-bath interaction is char- 
acterized by the spectral function 



g v (w) = 2tt^ \l >k 5{uj - uj k ). 



(4) 



In what follows we simplify the discussion and consider 
a two-level system (TLS) (n=0,l) for the local Hamilto- 
nian, 



H TLS 



APS copyright, to be printed in Phys. Rev. Lett. 



(*) 

yTLS 



B{t)a z /2- B(t) = [B + e(t)] 
F v a x . (5) 
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Here a x = |1)(0| + |0}{1|, a z = |1> (1| - 10)(0|, B = e{ 0) - 
Eq°\ and e(t) is the stochastic modulation of this energy 
gap. We sometimes refer to the truncated harmonic mode 
as a " spin" . In the weak coupling limit the TLS dynamics 
can be trivially carried back to the harmonic limit (n 
0. 1, ...oo), as we discuss below. 

Stochastically averaged master equation. Closed ki- 
netic equations for the \n) state population can be 
obtained in the weak system-bath coupling limit [l6l |. 
Briefly, starting with the quantum Liouville equation, 
this involves (i) studying the quantum dynamics within 
second order perturbation theory, (ii) going into the 
markovian limit, assuming short correlation time of bath 
fluctuations, and (ii) averaging over the stochastic pro- 
cess under the decoupling approximation, 



(j So P n ( x )) £ ~ (e^o^')dt'^ {Pn{x))e 



(G) 



Here P n (n = 0,1) denotes the TLS population, and the 
average is performed over energy fluctuations. The de- 
composition ([6]) relays on the separation of timescales: 
Energy fluctuations are assumed to be fast in compari- 
son to the characteristic subsystem relaxation time. The 
three assumptions (i)-(iii) result in the following equa- 
tions (f\(t)} e + <P (*))e = l, 

(A(*)>« =-(kLo + k^ )(Pi(t))e 

+ +*£ + l)<* , 0(t)>e, (7) 

with the activation and relaxation rates 



&o_^i = / dLjg u {oj)N u {oj)I(Bo—oj), 

J — oo 

/oo 
diog u (uj)[N u (oj) + l]I(B -oj). (8) 
-oo 

Here g„(w) is the phonon spectral function defined in ([J}, 
extended to negative values g v {io) = — <?„(— cj), N^(lo) = 
[e 13 "^ — I] -1 is the Bosc-Einstein occupation factor and 



1 r°° 

1(c) = — / dte iu)t {e i ^< t '^ dt ') 

27T J-oo 



(9) 



is the spectral lineshape of the Kubo oscillator. Assuming 
a Gaussian process, going to the fast modulation limit 
(e(ii)e(i2)) = 27<5(£i — tz), the noise spectral lineshape 
can be approximated by a Lorentzian function 



J( W ) = 



7A 
uj 2 + 7 2 



(10) 



of width 7. For 7 — > 0, I(u) = S(oj), and the field-free 
vibrational relaxation rates [17| are retrieved from Eq. 
©: The rates are evaluated at the bare TLS frequency, 
*JU = 9u(B )N u (B ), k^o = _gv{B )[N u {B ) + 1], 
obeying the detailed balance condition. In contrast, for 
a general driving this condition is violated, as the time 
dependent field drives the subsystem out of thermal equi- 
librium with the bath. 



The master equation ([7]) can be trivially generalized in 
the weak coupling regime to describe a harmonic local 
mode [Eqs. ([T])-([3])] under a stochastic field, 

(P„)e = + l)K^n+l + ^JUn-lK^e 

V 

V V 

(11) 

Note that for harmonic systems the microscopic rates 
are independent of the index n, thus we only identify 
activation and relaxation processes using ([8]). 

Heat current. We derive an analytical expression for 
the heat flux in the presence of stochastic energy mod- 
ulations. We begin with the TLS model ©, then ex- 
tend the results to the harmonic case. An expression 
for the energy flux operator, applicable in both station- 
ary and time dependent situations, was given in [l8l ]. 
j R = 1 [{Hj LS - H B ),V£ LS }. The current is defined 
positive when flowing from L to R, see Fig. [T] We sub- 
stitute Eqs. (J2J> and (O into this expression and obtain 

Mt) = -~ [<r*G R + B{t)a y F R ] , (12) 

with o y = -t|l)<0| +i|0)(l|, F v = Ek^M b l,k +Kk), 
G v = i^2 k K,k^k{b vk - b v . k ). The expectation value 

of the current is given by J = Tr{pJ}, where p is the 
total density matrix, and the trace is performed over 
the thermal bath and the subsystem degrees of free- 
dom. A further averaging over the TLS energy fluc- 
tuations is necessary. Since in steady-state the condi- 
tion (Tr{dV R /dt}} e = translates to (Tr{G R a x }) e = 
(Tr { B (t) a y F R }) e , the stationary noise-averaged current 
can be calculated by studying the first term in (TT2"j) . 
(Jr)* = -(Ttb{pi,qGr + po,iG R }) t . Here p i}j are the 
matrix elements of the density matrix, and Ttb denotes 
a trace over the bath states. Next, using a second order 
expansion of pi,o, adopting the same set of approxima- 
tions employed above for deriving ([7]), we acquire the 
steady-state heat current expression 



(jR)e = (Pl)ef^ ~ (Po)ef R 



(13) 



The noise-averaged population is obtained by solving ([7]) 
in steady-state, taking (P n )e = 0. The energy absorption 
and dissipation rates (Energy/time) are given by 

/oo 
dujug„(uj)I(B - u)N v (oj); 
-OO 

fUo = / dwug v (uj)I(B - uj)[N„(uj) + 1}. (14) 



Note that in the absence of fluctuations the rates sat- 
isfy the relation /^ n±1 = B k^ n±1 [T^], while in the 
presence of stochastic force we cannot trivially connect 
the energy transition rates (|14|) to the population relax- 
ation rates ([5]). In the weak coupling limit Eq. CG 
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can be easily generalized to the harmonic (Kubo) case, 

(J*)< = T,Zl n K P n)J^n-l - (Pn-l)ef*-^ n }, with 
the population obtained by solving (fTT|) in steady-state. 




FIG. 1: Scheme of the stochastic pump. A molecular element 
under an external perturbation (curly line) is coupled to two 
solids of different spectral properties. The arrows display the 
sign notation. 



Pumping mechanism. The pumping mechanism relays 
on two conditions: (i) The spectral function (of at least 
one bath) should strongly vary in the frequency window 
of the noise lineshape. (ii) The phonon reservoirs should 
have distinct spectral properties. We assume that lo r < 
Bq < u)l) lo v is the v contact cutoff vibrational frequency. 

Due to the fast stochastic modulation, the spacing 
B{t) can be abruptly reduced to values below u)r. Since 
these modulations are faster that the subsystem relax- 
ation time, the TLS effective temperature becomes signif- 
icantly low, even lower than Tr. The TLS then absorbs 
energy from both baths. In contrast, fluctuations that 
largely increase the gap lead to a high effective internal 
temperature, while practically disconnecting the R bath 
from the system. If the internal temperature is higher 
than Tl, energy will dissipate from the Kubo mode to 
the left bath. Thus, in this asymmetric setup net energy 
can be injected from the R bath to the L bath even for 
Tl > Tr. Note that this principle is independent of the 
specific form assumed for the noise spectral function, the 
molecular properties, and the details of the baths. 

This mechanism is reminiscent to the idea proposed in 
[lfj| . with two crucial distinctions: (i) Here random per- 
turbations catalyze (maybe even undesirably) heat flow, 
rather that carefully shaped pulses, (ii) The pump oper- 
ates in the fast modulation regime, rather than adiabat- 
ically/quasiadibatically, leading to the Carnot efficiency 
under proper conditions. 

An exactly solvable model. We employ next the Ein- 
stein model for the phonon spectral functions, g v {w) = 
£ v S{u> — lo v ), = 2nXl from (g]). Practically, this sce- 
nario may be realized by filtering the contact's phonon 
modes [19 ]. Assuming that I{Bq + lu u ) <C I {Bo — uj v ) 
[2C| . the transition rates and (fhT| reduce to 

fcuUi = Z»I(B - u u )N u {ui v ); 

kUo = CuI{B - u u ) [N u {u v ) + 1] ; 

/o_>i = oj v £, v I{B - u v )N v {w v ); 

f^o = w„£„J(B - w v ) [N v {u v ) + 1] , (15) 



and the energy current at the v contact becomes [l5| 

(J„) e = u„T [N l {uj l ) - N R (w R )\ (16) 



where T = lI^l^nL 



for the Hamiltonian dHJ) — ((3j) 

and T = ^ ^Hi+^U)] for thc TLS modcl ©' Here 
7^(1^) = £ U I{B — lo u )- In both cases the R contact is 
cooled down when Nr{ujr) > Nl{ojl), i-e. for (3l^l > 
PrijJr. Defining 6T = T L - T R > 0, 5lo = io L - lo r > 0, 
the pumping condition translates to 



ST Slv 

Tr UJr 



(17) 



For reversible systems the net current is zero, and an 
equality is attained, 5T/Tr — Suj/ujr. Note that for 
ST = a unidirectional current develops merely based on 
the asymmetry lol > ojr. We calculate next the machine 
efficiency. In steady-state the work over time done by the 
external force is given by the difference (W) e = (<7r)c — 



{W) t = {ur - w L )T \N l {uj l ) - N r {wr)\ . 



(18) 



The averaged cooling efficiency 77 = — {J R ) t / {W) t (the 
minus emerges due to the sign notation) then becomes 



_ u R T R _ 

f] — ^ < ^rp — Vmax, 



(19) 



where the inequality was obtained using (I17|) . Thus, 
the maximal efficiency of our cooling device is ex- 
actly given by the Carnot efficiency, reached when the 
device works reversibly. This result is independent 
of system properties (spin vs. harmonic mode) and 
the details of the noise lineshape. Finally, we ver- 
ify that the average entropy production per unit time 
is always positive for an irreversible device, (cr) e = 

rfft-zA [Nr{ur) — Nl(ul)]. 
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FIG. 2: Color-map of the heat flux at the R contact, revealing 
a pumping regime (central oval); ujl = 200, ujr=3, Tl = 
Tr = 25. For these parameters (Ji) e (not plotted) is negative. 

Numerical results. We exemplify the operation of the 
heat pump by simulating Eq. (|13|) with the stationary 
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FIG. 3: A color-map of the heat flux at the R contact where 
(Jr) € is negative for 7 < 8 and T L - T R < 4. lo l = 200, 
lor=S, Bo = 15, T R = 25. 

population determined by solving ([7]) in steady-state. 
We adopt the Lorentzian lineshape (fl~0|) . and assume 
that both solids are characterized by the ohmic spec- 
tral function g v (ui) = A u ue~ u / U ". We take uj l — 100 
and uir ~ Bo ~ 10, Al^r = 1. The calculations are 
performed in dimensionless units. If the energy parame- 
ters are given in meV, the current divided by the factor 
h = 4.14 (meV ps) yields fluxes in units of meV/ps. 

First we assume that Tl — Tr, and demonstrate a 
unidirectional heat flow from R to L. Fig. [2] presents 
a color-map of the heat current at the R side. The flux 
at the central oval (7 > 2; 12 < Bq < 20) is negative, 
i.e. the R bath is cooled down. We can explain this 
behavior as follows. For the ohmic function with u)r = 3 
the R bath is practically disconnected from the TLS at 
energies of the order B(t) > 15. Therefore, for Bq ~ 15 
fluctuations that reduce the gap lead to an overlap with 
the R phonon spectral function, and the TLS can accept 
energy from this contact. Fluctuations that significantly 
increase the gap (B(t) > 15) imply a complete isolating 



of the R contact, as its spectral function does not overlap 
with the molecular vibrational frequency. Energy is then 
solely exchanged between the (hot) TLS and the L solid, 
resulting in a net heat flow. 

We further demonstrate pumping of heat against a 
temperature gradient. Fig. [3] shows that (Jr)^ is neg- 
ative for 7 < 8 and T L - T R < 4. We verified that the 
total entropy of the system is increasing. Note that the 
pumped current is rather small, yet the object of this cal- 
culation is to demonstrate the ubiquity of the pumping 
effect, rather than to develop optimization schemes. 

In summary, nano-objects may pump heat effectively 
by exploiting random fluctuations, given that the (dis- 
tinct) solids are characterized by frequency dependent 
spectral functions. Our work reverses standard percep- 
tions of microscopic heat pumps: (i) Directing heat at 
ST = 0, and pumping heat at finite ST are facile pro- 
cesses, taking place in response to a general noise line- 
shape, (ii) Instead of carefully shaping the external 
pulse operating on the subsystem, simple bath "shaping" 
(mode filtering) may lead to maximal performance. 

Temperature reduction of nano-level objects is impor- 
tant, e.g. for cooling electronics devices, and for con- 
trolling chemical reactions and molecular dynamics. The 
pumping mechanism described here is pertinent to other 
systems, as one could realize an analogous exciton heat 
pump [2l| | . The effect might be observed in various sys- 
tems suffering random noise, for example, in nanome- 
chanical resonators. Considering a double clamped nano 
beam, the natural frequency could be mechanically mod- 
ulated due to the beam's mass variation resulting from 
adsorption-desorption processes 22]. The resonator fre- 
quency could be also electrostatically tuned by (stochas- 
tically) varying the bias applied to a gate electrode [23j . 
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